Beyond width and density: stable carbon and oxygen isotopes in cork-rings provide insights of physiological responses to water stress in Quercus suber L.

As climate change increasingly affects forest ecosystems, detailed understanding of major effects is important to anticipate their consequences under future climate scenarios. The Mediterranean region is a prominent climate change hotspot, and evergreen cork oak (Quercus suber L.) woodlands are particularly climatically sensitive due to cork (bark) harvesting. Cork oak's drought avoidance strategy is well-known and includes structural and physiological adaptations that maximise soil water uptake and transport and limit water use, potentially leading to reduced stem and cork growth. Trees' responses to cope with water-limited conditions have been extensively described based on cork-rings width and, more recently, on cork-rings density, in dendroecological studies. However, so far, tree functional attributes and physiological strategies, namely photosynthetic metabolism adjustments affecting cork formation, have never been addressed and/or integrated on these previous cork-rings-based studies. In this study, we address the relation between carbon and oxygen stable isotopes of cork rings and precipitation and temperature, in two distinct locations of southwestern Portugal-the (wetter) Tagus basin peneplain and the (drier) Grândola mountains. We aimed at assessing whether the two climatic factors affect cork-ring isotopic composition under contrasting conditions of water availability, and, therefore, if carbon and oxygen signatures in cork can reflect tree functional (physiological and structural) responses to stressful conditions, which might be aggravated by climate change. Our results indicate differences between the study areas. At the drier site, the stronger statistically significant negative cork δ 13C correlations were found with mean temperature, whereas strong positive cork δ 18O correlations were fewer and found only with precipitation. Moreover, at the wetter site, cork rings are enriched in 18O and depleted in 13C, indicating, respectively, shallow groundwater as the water source for physiological processes related with biosynthesis of non-photosynthetic secondary tissues, such as suberin, and a weak stomatal regulation under high water availability, consistent with non-existent water availability constrains. In contrast, at the drier site, trees use water from deeper ground layers, depleted in 18O, and strongly regulate stomatal conductance under water stress, thus reducing photosynthetic carbon uptake and probably relying on stored carbon reserves for cork ring formation. These results suggest that although stable isotopes signatures in cork rings are not proxies for net growth, they may be (fairly) robust indicators of trees' physiological and structural adjustments to climate and environmental changes in Mediterranean environments.

Using gradient Forest to predict climate response and adaptation in Cork oak.

Climate change is impacting locally adapted species such as the keystone tree species cork oak (Quercus suber L.). Quantifying the importance of environmental variables in explaining the species distribution can help build resilient populations in restoration projects and design forest management strategies. Using landscape genomics, we investigated the population structure and ecological adaptation of this tree species across the Mediterranean Basin. We applied genotyping by sequencing and derived 2,583 single nucleotide polymorphism markers genotyped from 81 individuals across 17 sites in the studied region. We implemented an approach based on the nearest neighbour haplotype 'coancestry' and uncovered a weak population structure along an east-west climatic gradient across the Mediterranean region. We identified genomic regions potentially involved in local adaptation and predicted differences in the genetic composition across the landscape under current and future climates. Variants associated with temperature and precipitation variables were detected, and we applied a nonlinear multivariate association method, gradient forest, to project these gene-environment relationships across space. The model allowed the identification of geographic areas within the western Mediterranean region most sensitive to climate change: south-western Iberia and northern Morocco. Our findings provide a preliminary assessment towards a potential management strategy for the conservation of cork oak in the Mediterranean Basin.

Phellem versus xylem: genome-wide transcriptomic analysis reveals novel regulators of cork formation in cork oak.

Cork cambium (or phellogen) is a secondary meristem responsible for the formation of phelloderm and phellem/cork, which together compose the periderm. In Quercus suber L., the phellogen is active throughout the entire life of the tree, producing a continuous and renewable outer bark of cork. To identify specific candidate genes associated with cork cambium activity and phellem differentiation, we performed a comparative transcriptomic study of Q. suber secondary growth tissues (xylem and phellogen/phellem) using RNA-seq. The present work provides a high-resolution map of all the transcripts identified in the phellogen/phellem tissues. A total of 6013 differentially expressed genes were identified, with 2875 of the transcripts being specifically enriched during the cork formation process versus secondary xylem formation. Furthermore, cork samples originating from the original phellogen (`virgin' cork) and from a traumatic phellogen (`amadia' cork) were also compared. Our results point to a shortlist of potentially relevant candidate genes regulating phellogen activity and phellem differentiation, including novel genes involved in the suberization process, as well as genes associated to ethylene and jasmonate signaling and to meristem function. The future functional characterization of some of the identified candidate genes will help to elucidate the molecular mechanisms underlying cork cambium activity and phellem differentiation.

Differential DNA Methylation Patterns Are Related to Phellogen Origin and Quality of Quercus suber Cork.

DNA methylation is thought to influence Quercus suber cork quality, which is the main constraint for its economic valorisation. However, a deep knowledge of the cytosine methylation patterns disclosing the epigenetic variability of trees with different cork quality types is totally missing. This study investigates the hypothesis that variations in DNA methylation contribute to differences in cork cellular characteristics directly related to original or traumatic phellogen activity. We used MSAPs (Methylation Sensitive Amplified Polymorphism) to assess DNA methylation patterns of cork and leaf tissues of Q. suber adult trees growing in three cork oak stands. The relationship between the detected polymorphisms and the diversity of cork quality traits was explored by a marker-trait analysis focusing on the most relevant quality characteristics. Populations differed widely in cork quality, but only slightly in degree of epigenetic differentiation. Four MSAP markers (1.3% of the total) were significantly associated with the most noteworthy quality traits: wood inclusions (nails) and porosity. This evidence supports the potential role of cytosine methylation in the modulation of differential phellogen activity either involved in localized cell death or in pore production, resulting in different cork qualities. Although, the underlying basis of the methylation polymorphism of loci affecting cork quality traits remain unclear, the disclosure of markers statistically associated with cork quality strengthens the potential role of DNA methylation in the regulation of these traits, namely at the phellogen level.

Patterns and drivers of scattered tree loss in agricultural landscapes: orchard meadows in Germany (1968-2009).

Scattered trees support high levels of farmland biodiversity and ecosystem services in agricultural landscapes, but they are threatened by agricultural intensification, urbanization, and land abandonment. This study aimed to map and quantify the decline of orchard meadows (scattered fruit trees of high nature conservation value) for a region in Southwestern Germany for the 1968 2009 period and to identify the driving forces of this decline. We derived orchard meadow loss from 1968 and 2009 aerial images and used a boosted regression trees modelling framework to assess the relative importance of 18 environmental, demographic, and socio-economic variables to test five alternative hypothesis explaining orchard meadow loss. We found that orchard meadow loss occurred in flatter areas, in areas where smaller plot sizes and fragmented orchard meadows prevailed, and in areas near settlements and infrastructure. The analysis did not confirm that orchard meadow loss was higher in areas where agricultural intensification was stronger and in areas of lower implementation levels of conservation policies. Our results demonstrated that the influential drivers of orchard meadow loss were those that reduce economic profitability and increase opportunity costs for orchards, providing incentives for converting orchard meadows to other, more profitable land uses. These insights could be taken up by local- and regional-level conservation policies to identify the sites of persistent orchard meadows in agricultural landscapes that would be prioritized in conservation efforts.

Fragmentation patterns of evergreen oak woodlands in Southwestern Iberia: identifying key spatial indicators.

Mediterranean evergreen oak woodlands (composed of Quercus suber L. and Quercus rotundifolia Lam.) are becoming increasingly fragmented in the human-modified landscapes of Southwestern Portugal and Spain. Previous studies have largely neglected to assess the spatial changes of oak woodlands in relation to their surrounding landscape matrix, and to characterize and quantify woodland boundaries and edges. The present study aims to fill this gap by analyzing fragmentation patterns of oak woodlands over a 50-year period (1958-2007) in three landscapes. Using archived aerial imagery from 1958, 1995 and 2007, for two consecutive periods (1958-1995 and 1995-2007), we calculated a set of landscape metrics to compare woodland fragmentation over time. Our results indicated a continuous woodland fragmentation characterized by their edge dynamics. From 1958 to 2007, the replacement of open farmland by shrubland and by new afforestation areas in the oak woodland landscape surrounding matrix, led to the highest values for edge contrast length trends of 5.0 and 12.3, respectively. Linear discriminant analysis was performed to delineate fragmented woodland structures and identify metric variables that characterize woodland spatial configuration. The edge contrast length with open farmland showed a strong correlation with F1 (correlations ranging between 0.55 and 0.98) and may be used as a proxy for oak woodland mixedness in landscape matrix. The edge dynamics of oak woodlands may result in different patterns of oak recruitment and therefore, its study may be helpful in highlighting future baselines for the sustainable management of oak woodlands.